AD‘‘A103  265  NAVAL  RESEARCH  LAB  WASHINGTON  DC  F/6  20/4 

THE  GROWTH  OF  HOT  REDUCED  DENSITY  CHANNELS  IN  GASES  DUE  TO  TURB _ ETC{ 

AUG  61  M  RALEIGH 

unclassified  NRL-MR-4555  NL 


REPORT  DOCUMENTATION  PAGE 


READ  INSTRUCTIONS 
BEFORE  COMPLETING  FORM 


2  GOVT  ACCESSION  MO.|  1  ACCIViEnT'S  catalog  NUMaCA 


A  Title  f«n«  $u6m»i 


THE  GROWTH  OF  HOT  REDUCED  DENSITY  CHANNELS 
IN  GASES  DUE  TO  TURBULENCE  AND  HEAT 
CONDUCTION 


J  AUThOA/a^ 

Michael  Raleigh 


*  AEAFOAminG  OAGAni2ATiOn  name  ano  aooacss 

Naval  Research  Laboratory 
Washington,  DC  20375 


<  <  CONTAOLLING  OAEiCE  name  amO  aooacss 

Office  of  Naval  Research,  Arlington,  VA  22217 
Defense  Advanced  Research  Projects  Agency 
Arlington.  V 


M  MOniTOAInG  AGENCT  name  a  AOOAESSTIf  AlllArani  Irom  Cwilralllna  OlOcAJ 


Naval  Surface  Weapons  Center 

White  Oak,  MD  20919  ATTN;  Code  R401 


*•  TvPC  A  ACAiOO  COVCMCO 

Interim  report  on  continuing 
NRL  problem 

S  AEAroAMING  OAG  AEAOAT  NUMSEA 


I.  contract  or  GKANT  NUM0CRf'*; 


10  PROGRAM  ELeMCNT.  RROjCCT  TASK 
AREA  «  WORK  UNIT  NUMfCRS 

6 1 1 5  3N  JIO 1 1  -09^  1  ;47-0S7 1  -0-1 
I  and  61101E;0,-OR40AA  f 


12  AEAOAT  DATE 

August  24,  1981 


<>  NLNSCA  or  RACES 

19 


IS.  SECuAiTv  Class,  tn  a>i,  r»paMj 

UNCLASSIFIED 

ISA.  OECLASSIEICATION/OOANGAAOING 
SCnEOULC 


17  O<STRl0uTiON  STATCmCHT  ,'«<  fA«  •A«CP«ef  «rtr«p»d  fft  Block  20,  II  diticront  /rom  Rmpon) 


surrucmcntarv  notes 

Research  supported  by  the  Office  of  Naval  Research  and  by  the  Defense  Advanced  Research 
Prefects  Agency  (DoD)  ARPA  Order  No.  3718,  monitored  by  the  Naval  Surface  Weapons 
Center  under  Contract  N60921-81-WR-WO1 14. 


19.  KEY  WORDS  (Cccitinvc  on  rovoroo  *idm  if  nmccccmry  and  Identity  By  klock  nvmkot) 

Reduced  density  chaimel 

Turbulence 

Heat  conduction 


STRaCT  fConiltnta  on  roooeo' 


oary  md  Idontity  By  ktack  nyatbor} 


Turbulence  causes  hot,  reduced  density  channels  in  gases  to  dissipate  rapidly.  The  scaling  of 
channel  radius  with  time  is  predicted  using  the  avera^d  governing  equations  of  fluid  motion  and 
the  assumption  of  pseudo-self-similar  expansion.  Correlation  terms  in  the  governing  equations  are 
replaced  by  terms  proportioiul  to  an  eddy  diffushrity  and  including  the  average  fluid  properties. 
The  diffusive  growth  thus  predicted  is  compared  to  experimental  results  and  to  the  normal  heat 
conduction  limit. 
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THE  GROWTH  OF  HOT  REDUCED  DENSITY  CHANNELS 
IN  GASES  DUE  TO  TURBULENCE  AND  HEAT  CONDUCTION 


I.  INTRODUCTION 

A  hot  channel,  fornied  by  either  laser  or  ohmic  heating  of  a  gas,''^  expands  initially  due  to  over¬ 
pressure.  If  spatial  asymmetries  exist  in  the  initial  heating,  vorticity  is  created  and  degenerates  into  tur¬ 
bulence  throughout  the  channel^.  A  second  phase  of  growth  then  follows  the  atuinment  of  pressure 
equilibrium  as  outside  gas  mixes  into  the  hot  channel.  Viscous  effects  cause  the  turbulence  to  decay 
until  heat  conduction  becomes  significant.  As  the  gas  responds  to  this  heat  flow  a  final  growth  and 
filling-in  of  the  channel  occurs. 

We  predict  the  scaling  of  the  channel  radius  with  time  during  both  later  phases  of  growth.  To  do 
this  we  use  the  governing  equations,  the  assumption  of  pseudo-self-similar  expansion  and  in  the  first 
case  a  simple  'mixing  length"  model  of  the  turbulence. 

II.  GROWTH  DUE  TO  TURBULENCE 

In  turbulent  flow  a  pattern  of  chaotic  eddies  with  a  spectrum  of  sizes  is  superimposed  on  the 
mean  fluid  motion.  If  the  fluid  velocity  is  represented  as  a  mean  value,  v^,  plus  a  relative  fluctuation, 
v',,  such  that 

Vr  -  V,  +  v’,. 

then  the  rms  value  of  the  fluctuation  velocity  is  a  measure  of  the  intensity  of  the  turbulence 

V,  -  (v;')''^ 

Turbulence  embedded  in  a  uniformly  moving  fluid  will  rapidly  become  isotropic 

Vv  -  V,  -  V.  -  V. 
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We  shall  model  the  channel  as  entirely  filled  with  isotropic  turbulence  and  undergoing  a  cylindrically 
symmetric  average  expansion. 

Turbulence  transports  fluid  properties  by  small  scale  convection.  However,  this  effect  may  be 
modeled  as  a  large  scale  version  of  transport  at  the  molecular  level  and  treated  as  an  anomalous 
diffusion.  We  suppose  small  volumes  of  fluid  move  a  distance  /  (Prandtl  eddy  length)  before  breaking 
up  and  losing  their  identities.  The  product  of  this  "mean  free  path"  and  the  relative  velocity  is  known 
as  the  eddy  diffusivity  and  becomes  the  effective  mass  diffusivity,  kinematic  viscosity,  or  thermal 
diffusivity 

D  ”  v  a  “  /v. 

fSlightly  different  values  of  /  actually  apply  to  each  situation.) 

Experiments  show  isotropic  turbulence  has  particularly  simple  properties.^  If  a  uniform  flow 
passes  through  a  grid,  the  resulting  turbulence  becomes  isotropic  a  short  distance  downstream  from  the 
grid.  The  intensity  of  the  turbulence  decays  beyond  this  virtual  origin  according  to 

1  1 

where  the  second  step  follows  from  the  constant  mean  velocity.  The  eddy  length  however  grows 

/  —  vx  ~  v7. 

We  may  therefore  ascribe  constant  transport  properties  to  the  medium  as  long  as  the  turbulence  per¬ 
sists.  The  value  of  this  consunt  depends  however  on  the  initial  conditions  under  which  the  turbulence 
was  created,  in  this  case  on  the  grid  wire  size. 

We  may  take  spatial  averages  of  the  variables 

p-l/pr/K  (1) 

v-^jfvdK  (2) 

f-~/prdP  (3) 
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where  we  select  the  averaging  volume,  K  as  larger  than  the  scale  of  the  turbulence  but  smaller  than  the 
channel  scale.  The  form  of  the  governing  equations,^  when  expressed  in  terms  of  the  average  proper¬ 
ties,  is  found  by  application  of  Reynolds  procedure^.  For  example  we  substitute  in  the  continuity  equa¬ 
tion  as  follows 

■^(p  +  p’)  -f  V-[(v  v')  (p  -I-  p’)l  *  0  (4) 

ot 

and  then  average  the  entire  equation  resulting  in 

-t- V  (pv) -l- V  (^)  -  0.  (5) 

In  recognition  of  the  transport  effect  of  the  turbulence  (and  following  Boussinesq’s  momentum 
transport  treatment’)  we  make  the  following  replacement  for  the  correlation  term 

pV  “  -aVp  (6) 

with  the  result  that 


+  Vpv  —  aV^p  =  0. 
ot 

(7) 

The  conservation  of  energy  equation  may  be  treated  similarly  and  the  correlation  terms  replaced  by  a 

diffusive  one.  The  governing  equations  thus  become 

^  -1-  V  p  V  -aV^p  -  0  conservation  of  mass 

Of 

(8) 

PCy 

+  vV7j  +  paRT„V  v  -  V  -  jpcyaV  zj  -  0  conservation  of  energy 

(9) 

pf  -  p„T„  equation  of  state 

(10) 

This  set  of  equations  incorporates  a  number  of  simplifications.  The  Navier-Stokes  equation  (conserva¬ 
tion  of  momentum)  has  been  eliminated  by  assuming  constant  pressure.  A  symmetric  expansion 
creates  no  new  turbulence  and  the  energy  contained  in  the  existing  turbulence  is  small.  Therefore  an 
additional  equation  to  describe  the  time  variation  of  v  is  not  needed  nor  do  we  find  a  term  in  the 
energy  equation  to  show  the  heat  generated  by  the  decay  of  turbulence.  No  terms  proportional  to  the 
dynamic  viscosity  appear  because  the  expansion  is  symmetric.  The  heat  conduction,  k,  is  neglected 
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during  the  turbulent  phase  in  comparison  with  the  anomalous  conduction.  pc\a.  We  may  rewrite  these 
equations  in  forms  appropriate  to  a  cylindrically  symmetric  expansion  (hereafter  dropping  the  averaging 
bars). 


9£L 

3r 


+  —  -j-  (rpv)  - 
r  9r 


1  A 

r  dr  dr 


0 


(11) 


pc. 


df  dr 


.  or  *  S  /  > 

pT~  P„T„. 


rpL\  a 


11 

dr 


(12) 

(13) 


The  approximate  scaling  of  the  channel  expansion  with  time  may  be  inferred  from  these  equations 
by  assuming  the  radial  profiles  of  T-T^,  p—p,,,  and  v  each  expand  in  a  pseudo-seif-similar  fashion. 
(These  profiles  share  a  common  characteristic  radius,  S,  at  any  time  f  under  the  presumption  of  con¬ 
stant  pressure.)  This  assumption  allows  us  to  replace  the  derivatives  by  simple  proportionalities  when 
we  apply  these  equations  to  a  point  coinciding  with  the  "edge"  of  the  channel.  (See  the  Appendix)  The 
various  constants  of  proportionality  (a  through  e)  are  retained  when  we  make  these  substitutions  in 
order  that  additive  cancellations  among  terms  may  be  seen  to  be  exact.  The  equations  now  reduce  to  a 
set  of  scaling  relationships 


(p-pj  .  ^vp  (p-p») 

a - +  d-^  -  aa - 5 - -  0 

r  S  8^ 


pc, 


(T-rj  ,  .  (r-rj 

C - - -  -t-  dv - ;; - 


^  __v  (r-rj 

+  ep„RT„j  -  cpcy  a  - - -  0 


pT~p„T„. 


(14) 

(15) 

(16) 


The  three  scaling  relationships  would  appear  to  be  insufficient  to  solve  for  the  unknowns  T,  p,  v,  and  8. 
However  after  eliminating  v  and  p  as  we  solve  for  S(r),  (T—  r,)  divides  out  so  that  the  scaling  is 
independent  of  temperature  as  expected  for  a  self  similar  expansion.  The  solutions  that  may  be  found, 

8^-ar  (17) 

v-0.  (18) 

where  8  and  v  are  the  mean  radius  and  radial  velocity,  show  that  the  channel  "diffusesT  away  at  a  rate 

determined  by  the  anomalous  coefficient.  This  behavior  persists  until  the  turbulence  decays.  Such 
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behavior  may  be  seen  in  Fig.  1  where  the  variation  of  the  square  of  the  radius  with  time  for  an  ohmi- 
cally  heated  channel  in  air  is  shown. 

m.  GROWTH  DUE  TO  HEAT  CONDUCTION 


We  now  consider  the  expansion  of  the  channel  due  to  heat  conduction  after  the  turbulence  has 
decayed.  Dropping  the  anomalous  diffusion  terms  and  retaining  the  thermal  conductivity  the  scaling 
relationships  become 


(T-TJ  ,  ,  (T-T„) 
c - +  flV - - - 


For  the  moment  we  treat  k  as  constant.  Then  a  solution  may  be  found  for  S(/) 


(19) 

v  (T-T„)  „ 

ep„RT,^  -  ck—^  =  0 

(20) 

PoTo- 

(21) 

Po(c.  +  fR) 


(22) 


Thus  the  combined  heat  flow  plus  fluid  response  scales  in  the  same  manner  as  ordinary  heat  flow. 


We  must  now  account  for  the  variation  of  the  heat  conductivity,  k,  with  temperature.  For  a  gas  k 
is  independent  of  density  but  proportional*  to  vT .  The  channel  behavior  now  depends  on  initial  condi¬ 
tions  as  evidenced  by  the  failure  of  T  to  cancel  out  of  the  scaling  relationships  when  we  assume 

k  -  kVT.  (23) 

An  approximate  integral  of  the  energy  equation  exists  which  is  independent  of  the  previous 
results  and  which  relates  T  to  S  through  the  initial  conditions.  We  take  a  global  view  and  treat  T and  p 
as  representative  values  for  the  entire  channel.  The  initial  heat  input,  Q,  must  always  equal  the  work 
done  by  expansion,  W,  plus  the  extra  internal  energy  E.  The  work  done  is  the  pressure  times  the 
volume,  at  density  p^,  needed  to  accomodate  the  gas  expelled  to  beyond  radius  5. 
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H'  -  p„ 


(ttSV,!  ~  ffS’p) 


P„RT„  ffS’ 


1  - 


The  extra  internal  energy  is 


£  »=  pi\irb^  ( T—  T„)  =  p„  T„i\7r&‘ 


1  - 


Equating  the  sum  of  £  and  to  0  we  find 


c  2 
^min 

5^ 


where 


(24) 


(25) 


(26) 


ST  =  _ ^ _ 

T  ■ 
^Po  '  n^p 


ill) 


(As  r  —  oo,  8  —  S^in  which  is  the  radius  at  which  a  channel  of  zero  density  and  infinite  temperature 
would  contain  Q  J/m  and  still  be  at  pressure  p„.)  The  previous  scaling  relations,  including  the  behavior 
of  k  with  T  imply 

ifVI, 


Po^p 


Eliminating  T  between  Eqs.  (26)  and  (28)  we  find  the  scaling  is  now 


(28) 


8V8'  -  5 I 
Pof-p 

where  k„  -  k^/T^. 


Under  normal  circumstances  8  >  5„in  and  the  scaling  reduces  to 

8^-  f 

Po^p 


(29) 


(30) 


This  limiting  slope  would  appear  as  horizontal  on  Fig.  1*. 


*A  ^meuk  factor  which  dependa  on  the  particular  jelf-fimilar  shape  and  choice  of  an  'ed(e'  has  to  this  point  been  assimilated 
into  a  For  a  Gaussian  proflie  with  the  I/e  point  as  the  edie.  equations  (17)  and  (30)  would  be  of  the  form  8^  -  Ax. 
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IV.  APPLICATION  TO  OHMICALLY  HEATED  CHANNELS  IN  THE  ATMOSPHERE 

Figure  1  displays  data  (8./)  for  an  ohmically  heated  channel  that  has  been  used  at  NRL  to  study 
the  propagation  of  relativistic  electron  beams  in  pre-formed  channels’.  These  channels  are  formed  by 
guiding  an  electric  discharge  from  a  Marx  generator  (K~250kv.  I  —  \QkA)  with  a  laser-induced,  air 
breakdown  and  have  been  produced  with  lengths  up  to  ~2m.  The  ohmic  deposition  is  ~3  J/cm  and 
the  channels  expand  to  a  ~1  cm  radius  in  —IQtis  before  they  reach  pressure  equilibrium  and  stabilize 
in  size.  After  —  150/z5  turbulence  becomes  evident  in  Schlieren  photographs  of  the  channels  and  they 
simultaneously  begin  the  growth  evidenced  in  Fig.  I.  The  kinetic  energy  in  the  turbulence,  ~10“^ 
J/cm,  is  a  small  fraction  of  the  internal  energy  of  the  channel.  It  nevertheless  strongly  influences  the 
channel  behavior.  Rapid  channel  growth  continues  for  —2ms  with  an  effective  dilfusivity  of 
—  500  cm^/ sec  which  is  ~  2000  times  the  thermal  diffusivity.* 

The  temperature  of  the  channel  strongly  affects  its  electrical  conductivity  and  thereby  the  propaga¬ 
tion  of  an  electron  beam.  In  these  experiments  the  electron  beam  was  expelled  from  the  channel  at 
t—  100/xs  but  showed  slightly  enhanced  propagation  in  a  channel  at  /  —  SOOfis  which  demonstrates  the 
importance  of  this  anomalous  cooling  mechanism. 

V.  CONCLLSIONS 

The  treatment  of  turbulence  as  a  diffusive  phenomena  is  often  used  in  engineering  analyses  to 
yield  useful  semi-empirical  results.  Reasonable  experimental  evidence  exists  for  applying  a  similar 
approach  to  the  cooling  of  a  turbulent,  reduced  density  channel.  The  gross  behavior  of  the  channel  is 
thus  "explained".  It  is  implied  for  example  that  the  radial  density  profile  of  a  Gaussian  channel  with  ini¬ 
tial  central  density,  p,,  and  initial  characteristic  radius,  6,.  evolves  in  time  according  to 

(8  1^  N' 

(p-po)- (p,-po)|y|  (31) 

_  7 

*ln  comparing  the  turbulent  expansion  to  true  heat  flow  a  Gaussian  channel  profile  is  assumed  and  allowance  is  made  for  the  ac¬ 
companying  factor  of  4  in  Equation  <1"). 


M.  RALEIGH 


where 

8^  =  8,'  +  4ar  (32) 

The  value  of  the  eddy  dilTusivity,  a,  remains  however  an  empirical  number.  The  treatment  presented 

here  assumes  a  constant  specific  heat.  c\  .  Such  an  assumption  is  valid  for  channels  in  air  provided 

T  <  3000  TC. 
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VTI.APPENDIX:  THE  RADIAL  AMD  TEMPORAL  GRADIENTS  OF  SELF  SIMILAR 
SOLUTIONS 

The  radial  profile  of  any  variable,  K  (eg.  T-  r„,  p-p„,  v,  pv)  is  here  assumed  to  expand  self- 
si  mi  larly 

.33) 

Where  simple  scaling  with  time  is  assumed  for  both  the  overall  amplitude  /  and  the  scale  length  8  such 
that 


f(i)~Pr".  (34) 

S(r)  -  yr.  (35) 

and  the  shape  is  always  that  of  g(xh  The  "edge”.  Xe^  is  a  fixed  value  defined  by 


where  e  is  a  fixed  fraction. 


yix?)  “  €g(0) 


(36) 


Consider  first  the  radial  derivative  of  K  at  a  point  r.  at  the  moment,  r.  when  the  "edge"  passes  by 

Br 


8 


hr-'.- 


^  Bx  Br 


hr"'.- 


/(r)|^ 


k.- 


1 

8(/) 


(37) 
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However  the  slope  of  at  the  edge  is  proportional  tog(x,  Vx.' 


(38) 


Thus 


dr 


8(t)  X,.  Sir) 

i.e.,  the  spatial  derivative  is  proportional  to  the  edge  value  of  •''divided  by  the  radius  ot  the  edge. 


(39) 


Consider  now  the  temporal  derivative  of  f 


d£ 

dr 


= 

dt 


+  g 


M 

di 


r  ■ 


Starting  with  the  first  term 


'dr 


fis,  d^ 

^dx  dr 


'  X.- 


Xr  or 


where  we  have  used  Eq.  (38). 

However  from  Eq.  (35) 

-7? 

Thus 

^X 

— m/- 

1 

dr 

7-'.“ 

If-. 

-mxv 

r 


The  first  term  thus  becomes 


(40) 


(41) 


(42) 


(43) 


fdg 

■'dr 


.  y(r.t) 
-  - 

ix  -  «.■  f 


(44) 


We  now  consider  the  second  term.  In  light  of  Eq.  (34) 


g^ 

^dr 


7  ■  r. 


g^nr 


/r- ! 


,  .  /■(/)  y(r.r) 

ng^Xe^  -  rr - - — 


(45) 


Combining  the  results  for  the  first  and  second  terms  we  find 
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ii! 

Br 


r  ■ 


,  yir.t)  .  I ■(/-,/) 
—  mi - +  n - 


yirj) 

l 


(46) 


so  that  the  temporal  derivative  is  proportional  to  the  edge  value  of  f  divided  by  the  elapsed  time  r. 


The  following  proportionalities  are  used  in  the  main  text 


dr 


a(p-po)  <p-p„) 

Br  / 

(4/1 

r  dr  0 

(48) 

MT  (T  -  Tn) 

Br  t 

(49) 

HT  .  (T  -  Tn) 

Br  ^  B 

(50) 

13.,  V 

—  — (rv)  =  e  — . 
r  Br  6 

(51) 

The  coefficients  for  the  second  derivatives  are  determined  by  the  requirement  that  the  equations 
reduce  to  the  proper  mass  and  thermal  diffusion  results  in  the  absence  of  the  terms  reflecting  fluid 
motion 


a 


ii. 

r  Br 


=  aa 


<P  -  Pq) 
6^ 


r  Br 


rpc\a 


Br 


cpL\a 


(/•-  To) 
6^ 


The  various  coefficients  are  interrelated  through  the  equation  of  state, 
equation  we  And 


(52) 

(53) 

Taking  derivatives  of  this 


Bp  _  ~Po^i)  BT 
Br  ~  p  Br 

Bp  ^  ~Po^o  BT 
Br  ”  p  Br' 

Dividing  Eq.  (54)  by  Eq.  (55) 

ie. 

Br  Br 

Bp  BT 

Br  Br 

but  in  terms  of  our  proportionalities 


(54) 

(55) 


(56) 


10 


NRL  MEMORANDUM  REPORT  4555 


SO  we  find 


a 


b 


-  Po> 
t 

(p  -  po) 
8 


c 


<  r  -  T-,,) 

t 


,  <  r  -  r„) 
d - ; - 


a</  “  be 

which  we  employ  in  reducing  the  scaling  relationships. 
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